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V!. DIGITAL COMPUTER SIMULATION OF THE DIGITAL IPS SYSTEM 
WITH WIRE-CABLE AND FLEX-PIVOT NONLINEARITY 

In this chapter the digital IPS with wire-cable and flex-pivot non- 
linearity is simulated on the digital computer. The main objectives of 
these simulations are to determine the effects of varying the sampling 
period T on the system stability, and to determine a range of "optimal" 
values of the parameters of the digital controller, K p , K. and 

The listing of the computer program is shown in Table 6-1. The 
parameters of the Dahl model for the flex-pivot nonlinearity are: 

Y = 32kbk 


T ppo = 2.25 x 10 N-M 

For the wire-cable nonlinearity, two ranges of values were used 

! wt : 

k WT - 100 ‘Vr ' °' 25 


for K^ t and H UT : 


h mt” 1 


H W t - 0-0' 


The nominal values of the digital controller parameters are. 


K = 70000 

P 

K. = 1,1 x 10 


K f = 19700 


6.1 Nominal Controller Parameters With Variable Sampling Periods 

The first set of simulations involves the study of the effects of 

various sampling rates. Figures 6-1 through 6-5 show the time responses 

of 0. due to an initial condition of 0.(0) = lo" 7 , for the nominal values 
i 

of K , K. and K.; K^ T = 0.25 and - 0.01, and T vanes from 0.005 sec 
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to 0.1 sec. As shown in Fig. 6-1. when T - 0.005 sec. the sailing period 
ls sufficiently small that the response of the digital system approaches 
that of the continuous-data system. As the sampling period increases, the 
decay time of the oscillation increases, and the system is marginally 
stable when T - 0.075 sec (Fig. 6-6). Figure 6-5 shows that when T >0.1 
sec. the response of a, diverges when time increases, and the system is 

unstable. 

6.2 Effects of Varying and H WJ 

cmjHv of the continuous-data IPS 
It has been shown earlier in the stu y 

that the values of and H,, do not affect the system performance 

substantially. Figures 6-6 and 6-7 are reruns of the cases in Figs. 

6-6 and 6-5, respectively, except with > 100 and V - 1 . These 

figures show that varying the values of and H UT by the wide margin 

did not change the responses of the system in any significant way. 

6.3 Effects of Varying K p 

In this section, simulations were conducted with 

K = Nominal value 

r 

K . = Nomi nal val ue 


T =0.1 sec 

and K is varied between 0.25 times to 5 times its nominal value. 

Figures 6-3 through 6-13 show the simulation results for 0., due to an 
initial condition of 0,(0) - lo' 7 . Figure 6-8 shows that when K p is 1/6 
its nominal value, the system is unstable. As K p is increased, the system 
becomes stable. However, when K p is 5 times its nominal value, the 
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Figure 6-3. 
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Figure 6-9 
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Figure 6-10 
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oscillations will take a long time to decay. These results show that 
there is an "optimal" value for K p from the standpoint of less oscillations 
and fast decay in the time response. Figure 6-11 shows that when K p 
equals to 2.35 times its nominal value the response seems to be the best 
under the simulated conditions. 

6.4 Effects of Varying K r 

Once a seemingly optimal value of K p has been determined, it is 
necessary to investigate the effects of varying the other parameters of 
the controller. The time responses shown in Figs. 6-14 through 6- 1 8 are fo 

K. = Nominal value 

K = 2.35 x Nominal value 
P 

T =0.1 sec 

*WT- ' 00 
h wt = 1 

and K varies between 0.5 times to 3-5 times its nominal value, 
r 

The time responses in these five figures show that as the value 

of K increases, the oscillations are reduced. However, for large values 
r 

of K , the response is slower. It may seem that an appropriate lange 
r 

of K is between two times nominal to three times nominal. The steady- 
r 

state error is also small with these values of K^. 

6.5 Effects of Varying K. 

Having selected the appropriate values of K p and K r in the initial 
round of t r i al -and-er ror design, the last parameter of the control let to 
be investigated is K. . In this section the effects of varying K. are 
investigated. Since the "optimal" value of K r seems to lie between two 
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times nominal and three times nominal, two sets of simulations were run, 

one with K = 2 * nominal and the other with K = 3 x nominal. Figures 
r r 

6-19 through 6-22 illustrate the time responses of 0. with 

K = 2.35 x nominal 
P 

K = 2 x nominal 
r 

T =0.1 sec 
Kyi - '00 

and K. varies between 0.5 times nominal to 2 times nominal. Figures 6-23 
through 6-28 show the time responses for K r = 3 x nominal and K. varying 
between 0.125 times nominal to 2 times nominal, with all other parameters 
the same as above. 

These simulation results indicate that for the sampling period of 
0.1 sec the system is not too sensitive to the variation of K.. In general, 
as K. increases the response becomes more oscillatory. Although lower 
values of K. give less overshoot, the steady-state error is actually 
larger. Therefore, a compromise has to be made between overshoot and 
steady-state error when making the selection on K.. Comparing the responses 
in Figs. 6-23 through 6-27, it seems that K. = 1.25 times nominal may be 
chosen from the standpoint of all-around performance. 

Figure 6-28 shows the time response of 0 . for the following param- 
eters : 
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The response Is extended for the time period up to 9.7 seconds to show 
that in the steady-state practica.ly no oscillations take place. 

t.6 Effects of Varying The Sampling Period 

Thus far ail the simulations have been carried out with a sampling 

period of 0.1 sec. Although the original system with the nom.nal values 

i i T — n 1 sec Fig. 6-28 shows that the 

for K , K. and K r is unstable w.th T - 0.1 sec, 9 

P ' • ri H hf> Stabilized and yield satisfactory responses 

digital IPS can indeed be stabilize 

by fine tuning the controller parameters. 

It is of interest to investigate if the digital IPS can be stabi.ized 

wlt h the given controller configuration with higher samp.ing periods. 

Figures 6-29 and 6-30 i.lustrate the 0, responses for T - 0.2 sec and 

0.3 sec. respectively. The other relevant system parameters are: 

K = 2.35 times nominal 
P 

K =3 times nominal 
r 

K. = 0.5 times nominal 


i k K and K. were chosen so that the system 

Note that the values of K p , K r and , 

performs best with T - 0.1 sec. Figure 6-29 shows that when T - 0.2 sec 
the system is stabie. a, though the response is siightiy osciiiatory during 
the transient state. However, the system is unstable for T - 0.3 sec 

(Fig. 6-30) . 
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6.7 Design of the Controller Parameters for T = 0.3 sec (K r Varies) 

It is apparent that the controller parameters that are suitable 
for the sampling period of T = 0.1 sec may not be the best for higher 
sampling periods. The objective now is to find a set of controller 
parameters so that the system has a satisfactory response for T = 0.3 sec. 

We shall first fix the values of K p and K. and vary Figures 6-31 

through 6-35 show the responses of 0. when K r varies between 1.25 times 
to 2, times its nominal value, and 

K =2.35 times nominal 

P 

K. = 0.5 times nominal 

i 

*WT > '00 

h wt = 1 

T =0.3 sec 

It is important to note that while the system is stable over a 
rather wide range of values for K when T = 0.1 sec, when T — 0.3 sec, 
the system becomes very sensitive to the value of K r . Figures 6-31 and 
6-35 show that for both small (K r = 1.25 times nominal) and large (2 times 
nominal) values of « r the response is very oscillatory. The response in 
Fig. 6-35 is actually unstable. Figures 6-32 and 6-33 show that K r = 1.5 * 
nominal or 1.65 x nominal seem to be the best choices for the specific 
situation. The response in Fig. 6-3^ for = 1.75 x nominal becomes 
slightly erratic at the 5-second mark. The simulation results indicate 
that at least the digital IPS can be stabilized for T = 0.3 sec. 

6.8 Design of the Controller Parameters for T = 0.3 sec (K p Varies) 

Since smaller values of K. would result in slower response, it was 
determined to vary only the value of K p in order to improve on the system 
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response for T - 0.3 sec. Figures 6-36 through 6-39 show the 0, responses 
for the following system parameters: 


K. = 

i 

K 

r 

= 

H WT = 


0.5 times nominal 
1 . 65 times nomi nal 
100 
1 

0.3 sec 


and K varies between 0.75 times nominal to 3 times nominal. These 
simulation results show that the K p - 1-25 * nominal case is the best 

among the cases studied. 

Figure 6-60 illustrates the 0. response of the above case for a longer 

. i L.. „ i v in ' + 2 x 10 seems 

time interval. An oscillation amplitude between 2 10 

to exist. This extremely small oscillation may be due to the digital 

simulation which has an integration interval of 0.02 sec. 

Figure 6-41 shows a case that the system response seems to be 

diverging after a sufficiently long period of time. 


6.9 Effects of Larger Sampling Periods 


Figures 6-42 and 6-43 illustrate the system responses when the 

. n a nd 0.35, respectively. The system 
sampling period is increased to 0.325 and u 

is ^stable under these conditions. When the sampling period is 0.5 sec, 
many simulation runs with various controller parameters did not yield 
a stable system. Figure 6-44 illustrates a typical unstable response for 

T = 0.5 sec. 

6.10 Conclusions 

Due to the nonlinear nature of the digital IPS system, the design 
carried out in the chapter is basically a trial-and-error effort. The 
simulation results show that the system can have fairly well-behaved 


I 


response for T 0.1 sec. For T - 0.3 sec, the system seems to perform 
well with the following set of controller parameters: 

Kp = 1-25 times nominal 
= 0.5 times nominal 
K. = 1.65 times nominal 

Although it may be possible to stabilize the system for sampling 
periods greater than 0.3 sec, fine tuning the controller parameters 

to achieve this is quite difficult, and the sensitivity of the system 
may be poor. 
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Figure 6-38. 
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